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Abstract 

Left-right symmetric extensions of the Minimal Supersymmetric Standard Model can explain 
neutrino data and have potentially interesting phenomenology beyond that found in minimal SUSY 
seesaw models. Here we study a SUSY model in which the left-right symmetry is broken by triplets 
at a high scale, but significantly below the GUT scale. Sparticle spectra in this model differ from the 
usual constrained MSSM expectations and these changes affect the relic abundance of the lightest 
neutralino. We discuss changes for the standard stau (and stop) co-annihilation, the Higgs funnel 
and the focus point regions. The model has potentially large lepton flavour violation in both, left 
and right, scalar leptons and thus allows, in principle, also for flavoured co-annihilation. We also 
discuss lepton flavour signals due to violating decays of the second lightest neutralino at the LHC, 
which can be as large as 20 fb"! at ^ = 14 TeV. 
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I. INTRODUCTION 



Left-right (LR) symmetric extensions of the MSSM ("Minimal Supersymmetric extension 
of the Standard Model") automatically contain the correct ingredients to explain the observed 
neutrino masses and mixings. The right-handed neutrino superfield, u'^, is necessarily part 
of the theory and breaking the LR symmetry by SU{2)ji triplets generates at the same 
time a Majorana mass term for the u'^ and thus a seesaw mechanism in contrast to 

type-I seesaw models where a Majorana mass term is merely added by hand. In LR models 
new gauge and Higgs fields appear with masses below the GUT scale, which change the 
running of all parameters under RGE evolution. In particular, the change in the running 
of the soft supersymmetry breaking masses can lead to potentially interesting effects in the 
phenomenology of SUSY models, even if the scale of LR breaking is above the energy range 
testable by accelerator experiments. 

Quite a large number of different LR models have been discussed in the literature. The 
original (non-supersymmetric) LR models jsj-Q break the SU{3)c x SU{2)l x S'f/(2)/j x 
U{1)b~l group by scalar doublets. Later it was realized that breaking the LR symmetry by 
(a pair of) triplets with B — L = 2 generates automatically also a seesaw mechanism j2|, |8|. 
Four triplets are needed jsl in the supersymmetric version of this LR model: A(l,3, 1,2) 
and A^(l, 1, 3, —2) due to the LR symmetry and A(l, 3, 1, —2) and A^fl, 1, 3, 2) for anomaly 
cancellation. Aulakh et al. 10|, lUl extended this minimal LR model □ introducing an addi- 
tional pair of triplets 3, 1, 0) and 1, 3, 0) with zero lepton number. The advantages 



of this setup, in the following called the fiLR model, are two-fold First, the LR break- 
ing minimum leaves R-parity unbroken already at tree level and, second, a non-trivial CKM 
matrix for quarks is generated easily and without resorting to flavour violating soft terms 
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Both, the minimal SUSY LR model Q as well as the fiLR model 10, 11| break the LR 
symmetry at an energy scale far above the range accessible to accelerator experiments. Only 
indirect tests of these models are therefore possible, all of which require some assumptions 
about the high scale boundary conditions for the soft SUSY breaking terms. Assuming 
CMSSM ("constrained" MSSM) [3l boundary conditions, two kind of indirect signals are 
possible, in principle: (i) Lepton flavor violating (LFV) decays and (ii) changes in the SUSY 
particle mass spectra. 

LFV decays are induced in supersymmetric models, even for strictly flavour blind bound- 
ary conditions, in the RGE running of the soft parameters, as has been shown for the case of 
type-I seesaw already in With LFV observed in neutrino oscillation experiments jl8| . 
one expects LFV to occur also in the charged lepton sector in practically all SUSY models. 
However, there is a qualitative difference between seesaw models and LR models: whereas 
in SUSY seesaw LFV is expected to occur dominantly in the left slepton sector 
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^ Several different realizations of LR models have been called minimal in the literature. These contain such 

Ii and Q,[li. While in the SUSY model of 



diverse variants as the models of 



introduced and B — L (and thus also R-parity) is broken by the vacuum expectations value of the i)'^, the 
non-SUSY models of (la. Il^ do not introduce any bi-doublet, but only a pair of left- and right- doublets. 
In this construction fermion masses have to arise ^rom non-renormalizable operators [isl . Il4| . 



an LR symmetry implies that left and right sleptons should have equal soft mass terms. 
Even with LR symmetry then broken at low energies, LFV in the right slepton sector can 



be sizable, as was shown in 
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Right slepton LFV could, in principle, be detected directly 



at accelerators or indirectly by measuring polarization in the decay /i — 67 (22[ 



Potentially measurable differences in SUSY mass spectra with respect to CMSSM ex- 



pectations can occur in seesaw type-II 
model 
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2J| and type-Ill 



24, 
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as well as in the f2LR 



21|. Changes with respect to CMSSM can best be understood analytically by form- 



ing certain "invariants", i.e. soft SUSY breaking mass parameter combinations which are to 
1-loop leading- log order constants over large ranges in CMSSM space j26|. Although there 
are quantitatively important 2-loop corrections to these invariants, as we will show below in 
the f2LR model the invariants have qualitatively different behavior from all seesaw models, 
due to the LR symmetry. This is similar to the situation discussed recently in 
context of different 5*0(10) based models. 
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Even slight changes in the SUSY spectra can lead to rather sizable changes in the calcu- 
lated relic densit y o f the lightest neutralino, f2^o/i^, as has been shown for seesaw type-II j28j 
and type-Ill 25|, |29| and also for an SO(IO) based model 30|. This fact can be easily under- 
stood taking into account that all solutions in CMSSM parameter space [3l|, which survive 
the latest WMAP constraints ji^l, require some special relations among SUSY masses in 
order to get a low enough Q^oh"^. For example, in the stau co-annihilation region VL^oh^ is 



essentially determined by Am 



m. 



>o, as long as Am < 10 GeV or so. Changing m^^ 



and/or m^o by only a few GeV will then change ^^o/?-^ by a large factor. Below we will 
discuss how the standard solutions to the DM problem are changed within the f2LR, with 
respect to the "pure" CMSSM allowed regions. Since the f2LR model has potentially large 
LFV in the right slepton sector, the DM constraint can also be fulfilled using the flavoured 



co-annihilation solution, recently discussed in 33|. We provide and discuss a few examples 



where flavoured co-annihilation can be realized in the parameter space of the VlLK model. 



Finally, LFV SUSY decays might be seen at the LHC. In 



21 we have shown that 



branching ratios xi ~^ hljXi^ with i ^ j are potentially large within the VtLK. Here we 
extend this work by calculating the SUSY production cross section for X2 from cascade 
decays, to estimate the number of events that could be found at the LHC. We have found 
a X Br{xl -> khx\) up to (10 - 20) fh. 



The rest of this paper is organized as follows. In the next section, we define the basics 
of the f2LR model. In section UTTl we discuss gauge coupling unification, generalities for the 
expected changes in the SUSY spectra with respect to CMSSM expectations and certain 
aspects of lepton flavour violation. In section HV] we discuss dark matter in the f2LR model, 
while section IVl discusses lepton flavour violation phenomenology for the LHC. We then close 
with a short summary. 
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Superfield SU{3)c SU{2)l SU{2)r U{1)b^l 
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TABLE I. Summary of the triplets of the ilLR model above the S'C/(2)r breaking scale. 



II. MODEL 



In this section we present the model originally defined in 10|, lllj, where a two step 



breaking of the LR symmetry was proposed in order to cure the potential prob 
to the conservation of R-parity at low energies 



34, 35 1 . For further details see 



ems related 
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1. Step 1: From GUT scale to SU{2)ji breaking scale 

Below the GUT scale@ the gauge group of the model is SU{3)c x SU{2)l x SU{2)r x 
U{1)b-l- In addition, parity is assumed to be conserved. Besides the quark and lepton 
superfields of the MSSM with the addition of (three) right-handed neutrino(s) u'^, some 
additional superfields are required to break the LR symmetry down to the standard model 
gauge group. First, two generations of $ superfields, bidoublets under SU{2)l x SU{2)ji, are 
introduced. They contain the standard Hd and Hu MSSM Higgs doublets. Note, that two 
copies are needed in order to generate a non-trivial CKM matrix at tree-level. Furthermore, 
triplets under (one of) the SU{2) gauge groups are added whose gauge quantum number are 
given in Table [H Note that the model contains B — L = ±2 and B — L = triplets. 

With these representations, the most general superpotential compatible with the symme- 
tries is 

>V = YqQ^Q^ + YlL^L" - + fLAL + fL^A^L" 
+ aAQA + a*A'Q^A^ + + a*f]"<l>$ 

+ MaAA + M^A^A^ + Mniin + M^n^il" . (1) 

Note that the superpotential in eq. ([1]) is invariant under the parity transformations 
Q O (g^)*, L ^ {L")*, $ o $t, A ^ (A=)*, A ^ (A^)*, n ^ (n^y. This discrete 
symmetry fixes, for example, the LPA'^L'^ coupling to be /*, the complex conjugate of the 
LAL coupling, thus reducing the number of free parameters of the model. 

Family and gauge indices have been omitted in eq. ([1]), more detailed expressions can 
be found in 10|. Moreover, the soft terms of the model can be found in 21|. The LR 



See subsection IIII Al for more details about gauge coupling unification and the GUT scale. 
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symmetry itself does not, of course, fix the values of the soft SUSY breaking terms. In our 
numerical evaluation we will resort to CMSSM-like boundary conditions at the GUT scale 
and obtain their values at the SUSY scale by means of the RGEs of the model. Finally, 
the superpotential couplings Yq and Yl are fixed by the low-scale standard model fermion 
masses and mixing angles. 

The breaking of the LR gauge group to the MSSM gauge group takes place in two steps: 
SU{2)r X U{1)b~l — ^ U{1)r X U{1)b~l — ^ U{1)y- In the first step the neutral component 
of the triplet Vf^ takes a VEV (vacuum expectation value): 

= ^ (2) 

which breaks SU{2)r. However, since /3ij(r2^°) = there is a U{1)r symmetry left over. 
Next, the group U{1)r x U{\)b-l is broken by 

(A^°) = ^ , (A^°) = ^ • (3) 



B-L 
2 ■ 



The remaining symmetry is now f/(l)y with hypercharge defined asY = I^n + 

The supersymmetric nature of the model is very relevant for the structure of the tadpole 
equations. Contrary to LR models without supersymmetry, the tadpole equations do not 
link fi^, and A'^ with their left-handed counterparts. Thus, the left-handed triplets can 



have vanishing VEVs [10| and the model leads only to a type-I seesaw. 

Neglecting the soft SUSY breaking terms and the electroweak symmetry breaking VEVs 
Vd and Vu (numerically irrelevant at this stage), and taking vbl = ^bl-i one finds the following 
solutions for the tadpole equations 

^« = ^ , VBL = -m^Mn)'/' . (4) 
ct a 

This implies that the hierarchy Vbl ^ requires Ma ^ Mq, as already discussed in flo| . 



2. Step 2: From SU{2)ji breaking scale to U{1)b~l breaking scale 

After the breaking of the LR gauge group to SU{3)c x SU{2)l x U{1)b x U{\)b~l some 
particles get large masses and decouple from the spectrum. In the case of the A triplets only 
the neutral components of the SU{2)b triplets, A^*^ and A^°, remain light. The charged 
components of the Vf^ triplet get masses of the order of Vji and thus only the SU{2)l triplet 
n and the neutral superfield ^l'^^ stay in the particle spectrum fill . [2l|. 

Similarly, the two bidoublet generations, $i and $2? get split into four SU{2)l doublets. 
Two of them remain light, being identified with the two Higgs doublets of the MSSM, 
whereas the other two get masses of order vr. As a result, the light Higgs doublets are 
admixtures of the corresponding flavour eigenstates, their couplings to quarks and leptons 
being combinations of the original Yukawa couplings, which has a strong impact on the low 
energy phenomenology. 
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The superpotential terms mixing the four SU{2)l doublets can be written as Wj\/ = 
{H{Y MhHI, where Hj^ = {H},H]^) and Hi = {Hl,H^) are the flavour eigenstates. In 
order to compute the resuhing couphngs for the hght Higgs doublets one must rotate the 
orig inal fields into their mass basis. Since Mh is not a symmetric matrix (unless = 0, see 



2l| ) one has to rotate independently and Hi, i.e. H^ = DH^, Hi = UH^. Here D and 
U are unitary matrices and H^ = (if^, if^) and HJ^ = {Hi, H!^) are the mass eigenstates, 
with masses m{Hjj^^) ~ (neglecting soft SUSY breaking terms and electroweak VEVs) and 
HHdJ ~ vr. 

The D and U rotation matrices are, in general, different. We can parametrize them as 

D - ( ^1 sin 6'i \ U - ( ^'^^ ^^"^ 1 (5) 

I — sin 9i cos 6*1 / ' I — sin 62 cos 62 J 

The angles 61 and 62 have a very strong impact on the low energy phenomenology. This 
can be easily understood from the matching conditions at the SU(2)r breaking scale. These 
include 

= — (sin^sld + sin^iFj , = ^(cos^2>d + cosO.Y^) , (6) 

■521 ^521 

Yl = — (sin e2Ye + sin OiY^) , Y^ = — (cos ^2^; + cos OiY^) . (7) 
■S21 ■S21 

where S21 = sin (^2 ~ (^i)- Eqs. Q and ([7j) show that in the limit 61 — t- 62 the high energy 
Yukawas Yq'^ and Y^'^ diverge. This implies that when 61 and 62 take similar values, 
lepton flavour violating effects at low energies, induced by RGE evolution by these Yukawa 
couplings, get larger. See section IIII CI for numerical results on this issue. 

Finally, neutrino masses are generated after the breaking of U{1)b-l through a type-I 
seesaw mechanism. Note that the superpotential in eq. ([T]) contains the term f*L'^A'^L'^ 
which, after U{1)b-l gets broken, leads to f^u'^u'^A'^^, with = — /*. Therefore, the 
matrix leads to Majorana masses for the right-handed neutrinos once gets a VEV. 
We define the seesaw scale as the lightest eigenvalue of the matrix Ms = JIvbl- 



III. GAUGE SECTOR, SPECTRUM AND LEPTON FLAVOUR VIOLATION 

In this section we discuss some important properties of the gauge sector of the model, 
the supersymmetric spectrum and the role of the bidoublet mixing angles for the size of the 
lepton flavour violating signatures. 



A. Gauge sector 

Despite the conceptual advantages of LR models, LR has fallen somewhat out of favor 
within the supersymmetric community. This is most likely due to the fact that within the 
MSSM gauge couphng unification is achieved automatically, if the scale of SUSY particles 
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jq12 jq13 jq14 jo15 jq16 jq12 jqB jqU jqI5 jgie 

VBL [GeV] VBL [GeV] 



FIG. 1. To the left, the gauge couphngs gsL = 92 (blue) and (red) at rfiQUT as a function of 
VBL- To the right, Aq as a function of vbl- In both figures vr is taken in the range [lO"*^^, 10"^^] 
GeV and rriGUT is defined as the scale at which qbl = 52- 



is of the order of 0{1) TeV or less. Additional intermediate scales tend to destroy this nice 
feature unless new particles and/or threshold effects are considered, see for example (36| . 

This issue is also present in the f2LR model. Let us define the scale ihcuT as the scale 
at which qbl = 92- In general, the strong gauge coupling does not unify with qbl = 92 
at ffiGUT and gauge coupling unification is not fully obtained. At 1-loop one finds: 

_ _ 87 55 

+ -^(29tBL + tR- msusv) (8) 
ozvr 

where tgusY = In {msusY/mz), Ibl = In {vBh/mz) and = In {vji/mz). Therefore, vbl is 
the most relevant scale for the determination of the difference between the gauge couplings 
at rfiGUT- This can be seen on the left side of Figure [H where the values of the three gauge 
couplings are shown as a function of vbl for vr € [10^"^, 10^®] GeV. For low values of vbl the 
running of qbl and g2 is too strong and their values at high energies clearly depart from the 
one obtained for g^. Furthermore, it is useful to quantify the deviation from unification by 
defining 

aBL{mGUT> + azymGUT) 

The right side of Figure [1] shows that for Vbl = 10^^ GeV the parameter Aq, can be as large 
as 0.45. Note, however, that this is not a principal problem. 

For instance, one possibility to recover gauge unification is the addition of new particles 
to the spectrum. As clearly seen in Figure [21 unification is not obtained due to a too fast 
running of 5^3. This can be fixed by adding new superfields charged under SU{2))c but singlet 



under the other gauge subgroups, as pointed out in |40|, |41[. Two examples are shown in 
Figure [3l where the addition of triplets of S'f/(3)c has been considered. To the left, one 
generation is added at msusY, whereas to the right five generations are added at vbl- In 
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l"gio (g?v) 

FIG. 2. 1-loop running of the gauge couplings for the choice of scales msuSY = 1 TeV, vbl = 10^^ 

2 

GeV and Vf> = 10^^ GeV. The dependence of , where = on the energy scale /x is shown. 
Different gauge couplings are represented in the different energy regimes. For /i G [tuzjVbl] one 
has 03 ^ (blue), aj^^ (green) and Oy^ (red). For ^ G [vbLj^b] one has 0:3"'^ (blue), (green), a]^^ 
(purple) and (red). For ^> one has a^^ (blue), a^^ = Oj^ = (green) and (red). 



both cases the new contributions to the running of (73 are sufficient to obtain gauge coupling 
unification. 



However, this picture might be a bit too simple: the authors of reference [40| pointed 
out that thresholds effects at the GUT scale can lead to important corrections to all gauge 
couplings. Note that these effects cannot be calculated unless a complete GUT model is 
specified and its high energy spectrum found. Since our model is motivated by an underlying 
50(10) gauge it is necessary to embedd also the A and f2 fields in complete representations. 
This will also lead most likely to a shift of g2 and g^L and the condition g2 = gsL using just 
the running values might lead a wrong prediction for the GUT scale. 

We fix this problem by using a fixed GUT scale mcuT = 2 ■ 10^^ GeV throughout the 
papeJl. In general, the three gauge couplings are different at this scale, but as long as vbl 
and vr are not much below mcuT the differences are not very large, e.g. of the order as in the 
usual MSSM. These differences would have important consequences in the gauge couplings 
and spoil known relations normally obtained by using unified gaugino masses. To take this 
into account we use gaugino masses at the chosen boundary scale of 2 ■ 10^® GeV which will 
unify at the correct GUT scale as explained in sec. IIIIBI Another advantage of this approach 
is that our GUT scale is no longer sensitive to the high energy VEVs vbl and vr. This has 
an important consequence: if we would have applied naively the universal SUSY breaking 



^ Note, that when using 2-loop RGEs one does in general not have strict unification 37H39| and one has 
to define a scale where the threshold corrections due to heavy degree of freedom have to be calculated, 
which is perfectly consistent with out approach. 
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2 4 6 8 10 12 14 16 2 4 6 8 10 12 14 16 



logio (gw) logio (g?v) 

FIG. 3. 1-loop running of the gauge couplings for the choice of scales msusY = 1 TeV, vbl = 10^^ 
GeV and Vfi = 10^^ GeV. Contrary to Figure [21 gauge coupling unification is obtained thanks 
to additional colored superfields. In the left panel, one triplet under SU{3)c, singlet under the 
other gauge subgroups, is added at msuSY, whereas in the right panel five generations of the same 
superfield are added at vbl- See Figure [2] for the color code. 



boundary conditions at mcuT, a large theoretical uncertainty in the determination of the 
low energy parameters would be present in our numerical conclusions. 

Let us now discuss another important theoretical issue: U{1) mixing. In the numerical 
implementation of the model we have chosen to break the left-right symmetric gauge sector 
down to the SM gauge group in one step in contrast to the heavy fields which we integrate 
out at two different scales. It would have been also possible to assume a gauge symmetry 
breaking in two steps, i.e. SU{2)ji x U{1)b-l — ^ U{1)r x U{1)b-l -> U{1)y- However, 

kinetic mixing. 



42 



Even if we 



the co-existence of two Abelian gauge groups introduces the possibility o 
because a term of the form nFj^^^F^'^'^^ is allowed by gauge invariance 
assume that this term were absent at the SU{2)ji breaking scale it would be introduced by 
RGE running already at 1-loop level because the two f/(l)'s are not orthogonal for the given 
particle content. This might be indeed relevant because it has been shown that the effect of 
kinetic mixing can be sizable even if the scale where the f/(l)'s coexist is rather short 
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Using the procedure presented in 4J| we have checked that the use of a two step breaking 
and running the parameters in the U{1)r x U{1)b-l basis is equivalent to the one step 
breaking and running the parameters under the resulting f/(l)y. Here, we had to use 1-loop 
boundary conditions and threshold effects to compensate differences in the 2-loop running 
as we use 2-loop RGBs. Nevertheless, although kinetic mixing is conceptually interesting, 
in the model under consideration its effects on the phenomenology are minor. 
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B. Low energy spectrum 



The introduction of additional superfields with masses below the GUT scale changes the 
RGEs not only for the gauge couplings but also for all MSSM soft terms. This in turn 
leads to changes in the electroweak scale supersymmetric mass spectra with respect to the 
standard CMSSM expectations. The CMSSM is defined by the following parameters: a 
common scalar mass mo, a common gaugino mass parameter M1/2 and common trilinear 
parameter Aq which are specified at the GUT scale. In addition tan /3 is specified at the 
electroweak scale and the sign of /i is fixed. As discussed in section 1111 Al we do not have 
exact unification of the gauge couplings but expect that threshold corrections can account 
for the difference. Also the gaugino mass parameters are subject to corrections of the same 
size. Therefore we define the boundary conditions for the gaugino mass parameters at the 
GUT scale mcuT as follows: 

M, = (10) 

OiBL 

For the calculation of the parameters at the SUSY scale we have used the complete 2-loop 
RGEs at all scales which have been derived with the Mathematica package SARAH 



45 



-47 



The numerical evaluation of the RGEs as well as the calculation of the loop corrected masses 

For more details 
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49| using the SPheno interface of SARAH 
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was done with SPheno 

about the numerical implementation see also 

An example spectrum is given in Figure H] where we show two SUSY (and Higgs spectra) 
for a specific choice of CMSSM parameters (mg = 600 GeV, M1/2 = 700 GeV, Aq = 0, 
tan/3 = 10 and fi > 0). Dashed lines are the mass spectra for the pure CMSSM case, 
while the full lines have been calculated for the QLK model with the specific choice of 
vr = vbl = 10^^ GeV. These parameters were chosen to lie outside the region already 
excluded by the SUSY searches at the LHC experiments 51-5^5 t>ut are otherwise arbitrary. 



These spectra should not be taken as predictions for the fiLR model, but serve only for the 
sake of discussing the main differences between this model and the pure CMSSM case. 
In general, with < < mcuTi the VtLK has a lighter spectrum than CMSSM for the 
same parameters. One can understand this semi-quantitatively with the help of the following 
considerations. 

Gaugino mass parameters run like gauge couplings do. Since in the f2LR model the 
running of the gauge couplings is changed with respect to the MSSM case, also gaugino 
masses are changed. Consider for example Mi. At 1-loop leading- log order we find 

Mi(m,^,,) = Mv2 ""^7^^^"i (11) 

= Mi/2 [Xi+X2(-3/i+/2)] 

where we have defined 

^ _5 3 ay(mg^gy) 15ay(mg^gy) / mcjjr \ ^ _ SayimsusY) 
^ 2 2 a2imsusY) 27r \msusY ) ' ^ 27r 
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FIG. 4. Example spectra comparing CMSSM with the modeL Parameters have been chosen as 
mo = 600 GeV, M1/2 = 700 GeV, Aq = 0, tan/3 = 10, /x > and vr = vbl = lO^'' GeV. FuU hnes: 
OLR, dashed hne pure CMSSM. Note that in this example one has rn{x2)cMSSM — 'm'{Xi)nLR by 
accident. 



and 



l,=in(^),k = ln(^). (13) 



Note that, due to their definition, /1.2 > 0, with /i = /2 = in the CMSSM Umit. One 
can easily check that Xi > and X2 > 0. For mcuT = 2 ■ 10^^ GeV and msusY = 1 TeV 
one obtains roughly Xi ^ 0.424 and X2 — 0.008. Therefore, Mi^rrisusY) decreases with 
li and increases with 12- In other words, one can decrease Mi{msusY) with respect to the 
CMSSM value by using a large li and /2 = and increase it only in the case I2 > 3/i. Similar 
equations hold for M2 and M3. However, at low energy ratios such as M1/M2 still follow 
the standard CMSSM expectations, only the relation to M1/2 is changed. 

Similarly, sfermion mass parameters can be written schematically as, ~ ml + CjMy^, 
where the coefficients are different for different sfermions and depend on vr and Vrl- 
Exact expressions are given in appendix \^ One can show that for vrl vr < mcuT the 
Cj are always smaller than in the CMSSM limit, explaining why also the sfermions in the 
f2LR model are lighter than in CMSSM, see Figure HI It is important to note, that gaugino 
masses change faster with vr than sfermion masses. We will come back to this point in the 
discussion about dark matter, see section IIVI 

Individual SUSY masses depend strongly on the initial values of mo and Mi/2- However, 
one can form four different combinations ("invariants") of soft SUSY breaking parameters, 
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10" 10 '5 10 



Vfi [GeV] 

FIG. 5. Invariant mass combinations QE/10 (pmple, dashed-dotted lines), LE (red, dashed lines), 
QU (dark blue, dotted lines) and DL/10 (light blue, solid lines) as functions of vr. vbl G [10^^, 10^^] 
GeV, upper lines correspond to vbl = 10^^ GeV and lower lines to vbl = 10^^ GeV. 



which we choose as 

LE = {ml - ml.) /Ml, 
QE^ {ml -ml) /Ml 

DL^{ml-ml)/Ml ^ ^ 

QU^ {ml -ml) /Ml 

where at the leading- log level m^ and M1/2 drop out. Formulas for an analytical 1-loop 
leading-log calculation of these four invariants are given in appendix 1X1 

Figure [5] shows the analytically calculated values of the invariants as a function of vr for 
two values oi vbl- The invariants depend more strongly on vr and only weakly on Vbl- As 
the figure shows, at leading order QE and DL show a very mild dependence on the scales vr 
and vbli while LE and QU decrease with decreasing vr. At the vr scale one has m\ = ml 
and ml = ml = ml due to parity conservation. This implies that the lower vr is, the 
smaller the difference between them at the SUSY scale is, which means that the LE and 
QU invariants can have values below the CMSSM prediction. We have checked that this 
effect is quite robust, and the theoretical uncertainties do not have much influence on it. 
In particular, these two invariants do not show any numerical dependence on mcuT, due to 
cancellations among left and right contributions in the running from the GUT scale to vr. 
This is an interesting result, since in CMSSM models with a type-II or type-III seesaw, all 
four invariants are always larger than their CMSSM limit for seesaw scales below the GUT 
scale 
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25 



|26|. This allows, at least in principle, to distinguish between models with a 
high scale LR group and CMSSM models with type-II or type-III seesaw as the explanation 
for the observed neutrino masses. 

Figure [5] shows the leading- log calculation. It is known in the case of seesaw based models, 



that there exist important 2- loop corrections and 1-loop threshold effects 23|, |25|. Similarly, 
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FIG. 6. Invariant mass combinations QE/10 and LE as functions of vbl (left) and vr (right). The 
color code is the same as in Figure [5j Solid lines correspond to the analytical approximation, dashed 
lines to 1-loop numerical results and dotted lines to 2-loop numerical results. Similar deviations 
are found for the other two invariant mass combinations. 

also in the QLK model important numerical corrections exist, as is shown in Figure El 
When the invariants are written including the effect of non-unification at the GUT scale, 
see appendix [XJ they reproduce qualitatively the numerical results. However, as shown in 
Figure [6] quantitatively important shifts are obtained going from leading-log to (numerically 
solved) full 1-loop calculation. Even going from 1-loop to 2-loop calculation, numerically 
important differences are found. 

To summarize, the invariants are good model discriminators in principle. Especially 
noteworthy, in the f2LR model LE and QU are expected to be below their CMSSM limit. 
However, to identify such spectrum distortions, once the SUSY mass spectrum is measured, 
will require highly accurate measurements and therefore measurements at an ILC. 



C. Lepton flavour violation and the role of 2 



Flavour violation in leptonic processes has attracted a lot of attention in the experimental 
community. Decays like /i — 67 have been searched for decades, without any positive result. 
Very recently, the MEG experiment 55| has published the results of the analysis of the data 
collected in 2009 and 2010 Q, setting the new bound Br(/i 67) < 2.4 ■ 10"^^ This 
impressive experimental limit strongly constrains models with extended lepton sectors, such 
as the QLR. model. 

The branching ratio for /j — ?■ Ijj can be written as (22I 



) Br(/, 



(15) 



The couplings Al and are generated at the 1-loop level. The relation between these 
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couplings and the slepton soft masses is given approximately by 



SUSY 



m 



(16) 



SUSY 



where msusY is a typical supersymmetric mass. In the derivation of this estimate it has 
been assumed that (a) chargino/neutralino masses are similar to slepton masses and (b) 
A-terms mixing left-right transitions are negligible. 

Note that, due to the negligible off-diagonal entries in m^c, a pure seesaw model predicts 



^ 0. On the contrary, in reference 2l| it was pointed out that a left-right symmetry 



at high energies induces non- negligible off-diagonal elements in m^c, giving additional con- 
tributions to LFV processes. In fact, taking into account the running from the GUT scale 
to the vr scale, the off-diagonal elements of the slepton soft masses at 1-loop order can be 



written in leading-log approximation as |21l . |5 



= ~ (Vf + E YfYp t) (3m? + Al) In {^-^^ (17) 

(ml.)y k = (mi.), |„. = - jij (3/'/ + |: rf 'tyf ) (iml + Al) In (18) 

which are of the same size in the CP conserving case. The A parameters also develop LFV 
off-diagonals in the running. We do not give the corresponding approximated equations 
because they do not lead to qualitatively new features. 

The angular distribution of the outgoing positron at, for example, the MEG experiment 



could be used to discriminate between left- and right-handed polarized states 
MEG is able to measure the positron polarization asymmetry, defined as 
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If 



A 



\A, 



\A 



R\ 



LR 



we will have an additional tool to distinguish from pure seesaw models, where A^r 
predicted. 



(19) 



1 is 



We extend here the discussion of ref. 
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addressing the influence of the high energy 



parameters on low energy observables. As explained in section [TTl the mixing angles of the 
bidoublets at the vr scale have a very strong impact on low energy LFV processes. As is 
shown in eqs. ([6]) and ([7]), this is due to their relation to the Yukawa parameters. If the 
bidoublet couplings are tuned such that 9i and 62 get very close values, the entries of the 
Yukawa matrix will become very large, which in turn will lead to large running effects for 
the soft squared masses of the sleptons. Therefore, we expect a strong correlation between 
S21 = sin(^2 — ^1) and the size of the LFV processes. 

The left side of Figure [7] shows such correlation for a particular but representative set of 
parameters. We note that a proper choice of 62 — Oi can enhance (or suppress) the /j — > /j7 
branching ratios by several orders of magnitude. Moreover, special values for 62 — 61 are 
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FIG. 7. To the left, Br(/i — t- 67) (blue), Br(T — t- /i7) (red) and Br(r — ?• 67) (green) as a function of 
I sin(02 — ^i)|- The horizontal line shows the curr ent MEG bound Br(/i 67) < 2.4-10-^2 ^j^g 
right, the positron polarization asymmetry Alr as a, function of | sin(02 — ^i)|- Both figures have 
been made fitting neutrino data with the flavour structure in and taking the parameter choice 
Ms = 10^2 Qev, VBL = VR = 10^5 GeV for the CMSSM point niQ = 120 GeV, M1/2 = 600 GeV, 
tan (3 = 10, Aq = and ^ > 0. 

found where the r decays have strong cancellations. Their origin can be traced back to the 
flavour structure of the LFV entries in and m^c. 

The RGE running from mcuT to vr introduces LFV entries in m| and m^c proportional 
to Ylk yj^'^^Yj^''^ ^ (the dagger is to be exchanged for the case of m^c, which makes no difference 
if CP is conserved). Expanding this expression in terms of Ye and Y,y, the Yukawa couplings 
at vr, one obtains 

E "^l'M'^ ' = ^ [^^^^ + ^'^^e + c.ayx + y.y})] (20) 
k ^21 

where C21 = cos(^2— ^i)- For the observable Br(/j — j- /j7) we have to compute (^^ y^'^)y^'^) t-j 
with i ^ j. The Y^Ye piece can be neglected, since Y^ is almost diagonal (we choose to work 
in the basis in which Y^. is diagonal at the electroweak scale). Furthermore, in the approxi- 
mation mi^ ^ mi. one can also neglect the Y^Y} piece, since it will be proportional to mi- . 
Therefore, using eqs. ffT^ and ( TTO]) . one finds 

Br(/, ^ /,7) cx ^\{YX + C21YXW (21) 

From eq. (12T|) one can see that as 9i and 62 get closer S21 — )■ and one gets an enhancement. 
The interesting new point is the interplay of the two terms as for specific values of 62 — 61 
cancellations occur. In case of of Br(/i — )■ ej) this hardly occurs as in general the //-Yukawa 
coupling is much smaller than the 1(1^,^)4^1 and a cancellation is not possible. On the contrary, 
the r- Yukawa coupling is comparable to Y^, in large parts of interesting parameter space, 
allowing for potential cancellations. 
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However, we want to stress that the position and the degree of these cancellations depend 
on the high energy VEVs vbl and vji and thus one cannot determine 62 — 61 uniquely. 
Moreover, this connection between leptons and bidoublet mixing angles is lost if one fits 
neutrino data using the superpotential coupling /, which gives rise to the Majorana mass of 



the right-handed neutrino. In this possibility, called '/ fit' in [2l|, there is no dependence 
on the bidoublet mixing angles, since / couples the leptons to the A triplets, but not to the 
bidoublets. For completeness we note that also the positron polarization asymmetry Alh 
depends somewhat on 62 — 9i as can be seen on the right side of Figure [71 

Finally, we have checked that this strong dependence on 61^2 is not present in the quark 
sector. Flavour violating processes, such as 6 — > S7, are very weakly affected by the bidoublet 
mixing angles and the phenomenology follows very closely the well-known results of the 
CMSSM. This is due to the fact that flavour violation in the (s)quark sector is dominated 
3y the CKM matrix. In other words, the f2LR model is a minimal flavour violating one 
3,01. 



IV. DARK MATTER 



Astrophysical observations and the data from WMAP 



32 



put on solid grounds the 

existence of non-baryonic dark matter in the universe. The PDG j62| quotes the value 
f^DM^^ = 0.110 ±0.006 at 1(7 CL. It is well-known that with this low value of ^dm^^ only a 
few, very specific regions in the parameter space of CMSSM can give the correct relic density 
for the neutralino 
region 



63| . These are known as (in no specific order) (i) the stau co-annihilation 



64| ; (ii) the stop co-annihilation region 65l-l67j: (iii) the focus point line 
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69 



and 



(iv) the Higgs funnel 

In any of the co-annihilation regions the relic density of the lightest neutralino, which in 
the above cases (i) and (ii) is mostly bino, is reduced with respect to naive expectations by 
a close mass degeneracy between the neutralino and the NLSP. Any changes in the SUSY 
spectra can then have a rather large impact on the calculated Q^oh"^. In the focus point 
region a comparatively small value of /i, with respect to the rest of the CMSSM parameter 
space, leads to an increased higgsino content in the lightest neutralino. This in turn leads to 
a larger coupling of the LSP to the Z^, thus reducing the relic density. The fourth allowed 
region, the Higgs funnel, appears for high values of tan /3 and for parameters where the CP- 
odd Higgs scalar A has a mass which is twice the LSP mass. In this case there is an s-channel 
resonant enhancement of the neutralino annihilation cross section and large values of tan 13 
are needed for this resonance to be effective, since the width of A is large at large tan (3. All 
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71 



of the regions discussed above have been extensively studied in the CMSSM 
and in some extensions of it, like CMSSM plus seesaw of either type-II or type-Ill 

Below we discuss the changes for these allowed regions for the f2LR model. We have 
used SPheno to compute the low energy spectrum and used it as an input for mi cr Omegas 



72| to obtain the value for Q^oh^. To include effects of flavour violation, we have created 



suitable model files for micrOmegas with SARAH. In the present setup the main effect on dark 
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FIG. 8. Stau co-annihilation. On the left: mass difference between the lightest neutralino and 
the lightest stau vs. vbl for three different values of vji = 10^^ (blue), 10^^ (red) and 10^^ GeV 
(green). On the right: resulting dark matter relic density Qh'^. The gray band shows the Qh'^ = 
[0.1018, 0.1228]. The other input parameters have been chosen as: uiq = 105 GeV, M1/2 = 600 GeV, 
tan/3 = 10, Ao = GeV. 

matter is via changes in spectrum (see section [III BP . For low Vbl ot Vr the dark matter 
allowed regions are reduced compared to the CMSSM. We have found this effect for the stau 
co-annihilation, stop co-annihilation, Higgs funnel and focus point regions and discuss these 
regions in turn. 

A. Non-flavoured dark matter allowed regions in riLR 

1. Stau CO' annihilation 

As discussed in section [III B| lowering the values of Vbl and leads in general to a lighter 
spectrum compared to the pure CMSSM case, for the same values of CMSSM parameters. 
In Figure [8] we show an example for the mass difference between the lightest neutralino and 
the lightest scalar tau for three values of Vr as a function of Vbl- Since the mass of the bino 
decreases faster than the mass of the stau, the mass difference increases and this in turn 
leads to a larger value of the neutralino density, as shown in the plot on the right, because 
co-annihilation becomes ineffective for = m^^ — m^o larger than a few GeV. Note that 
the CMSSM parameters in this example have been chosen to give approximately the correct 
relic density for vn ~ 10^^ GeV. 

Two examples of allowed stau co-annihilation regions are then shown in Figure [9l The 
figure shows the 3 a allowed regions for tan/3 = 10, and vbl = vr = 1.5 x 10^^ GeV (left) 
and Vbl = vr = 10^^ GeV (right). Lowering vr (and vbl) shifts the allowed region towards 
smaller values of tjiq. Depending on the values of tan/3 and Aq chosen, for values of roughly 
Vbl = Vr ^ 10^^ GeV the stau coannihilation line disappears completely. The observation 
of a SUSY particle spectrum consistent with the stau co-annihilation region could thus be 
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FIG. 9. Dark Matter allowed regions for tan/3 = 10, and vbl = vr = 1.5 x 10^^ GeV (left) and 
vbl = vr = 10^^ GeV (right). For this figure we have chosen a low seesaw scale, Ms = 10^^ GeV, 
which assures that the LFV constraints are respected. 



turned into a lower limit on Vbl and vr within CMSSM, at least in principle. 

However, as has been noted in section [IIIB[ the effect of vr on Mi is stronger and inverse 
to the one of Vbl- We have checked that it is possible to obtain a stau-coannihilation region 
slightly above the CMSSM expected region (for any fixed tan /3 and Ao) for values of vr 
close to the GUT scale and vbl low. Thus, the stau coannihilation can give a lower limit 
on vr only as a function of the (assumed) value of Vbl- 



2. Stop CO- annihilation 

With typical choices of CMSSM parameters the squarks are much heavier then the LSP. 
However, if the off-diagonal elements in the squark mass matrices are significant, the lighter 
mass eigenstates have their masses lowered and the lighter stop can be almost degenerate 
in mass with the lightest neutralino. This stop co-annihilation has been studied in the 
literature GS-G^I and happens for instance if the soft breaking parameter Aq has a large 
value We also explored this possibility in the present setup. Figure [10] shows an example. 
In this plot, the CMSSM parameters where chosen to be rriQ = 630 GeV, M1/2 = 400 GeV, 
tan/3 = 10, Aq = —2 TeV and /i > 0, which leads to a mass difference between stop and 
neutralino below 2 GeV in the pure CMSSM limit. Lowering vr or vbl below the GUT 
scale, this mass difference again increases, rendering co-annihilation ineffective. The lowest 
possible values for vbl and Vr in this scenario are usually found in the region vbl ~ vr, 
as the figure shows. We found no combination of parameters which improves the result 
shown in the figure by more than a few GeV. Our conclusion is therefore, that the stop co- 
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FIG. 10. m^o — m^^ contours in the vbl — vr plane. The CMSSM parameters take the following 
values: mo = 630 GeV, M1/2 = 400 GeV, tan/3 = 10, ^0 = -2 TeV and ^ > 0. 

annihilation region vanishes for low intermediate scales, similarly to the stau-coannihilation 
region. Observation of a SUSY spectrum consistent with stop co-annihilation could therefore 
be interpreted lower limit on vr and Vbl- 

3. Higgs funnel 

As in the previous two examples, we have found that also the Higgs funnel region depends 
rather strongly on the choice of and Vbl- As mentioned above, at large values of tan/3, 
typically tan /3 > 50 the width of the CP-odd Higgs boson A can be large enough so that 
there is a s-channel resonance for the neutralino pair annihilation in the region 2m^o ^ m^. 
An example for this Higgs funnel region is shown in Figure [TTl 

In the VtLK model the Higgs funnel region disappears when the intermediate scales are 
lowered below a certain limit. This is due to the existence of an upper limit on tan /3, stronger 
in the f2LR model than in CMSSM, caused by the requirement of (1) perturbativity of the 
bottom and tau Yukawa couplings, and (2) stability of the electroweak symmetry breaking 
minimum of the Higgs potential. In the CMSSM this limit is roughly tan /3 ~ 55 — 60. In 
the VLLK model this limit depends on Vb„ and becomes stronger when this high energy scale 
is lowered. This can be understood, in principle, from the fact that the lower vr is, the 
larger the ROE effects on the Yukawa couplings are. As shown in Figure [121 which should 
be taken as an illustrative example, for a given value of one can always find an upper 
limit on tan/3, which in the limit vr approaching mcuT = 2 ■ 10^^ GeV recovers the usual 
CMSSM limit on tan /3. For lower values of vr, therefore, there is an upper limit on the A 
width. The allowed Higgs funnel region becomes smaller until for a certain value of vr it 
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FIG. 11. The allowed region for ^y-oh'^ in the mo-Mi/2 plane for tan/3 = 53, vbl =vr = 10^^ GeV. 
For this figure we have chosen a low seesaw scale, Ms = 10^^ GeV, which assures that the LFV 
constraints are respected. 




vr [GeV] 

FIG. 12. Upper limit on tan /3 as a function of vr for vbl = 10^^ GeV. The figure was obtained for 
the choice of CMSSM parameters mo = M112 = 1-5 TeV, ^0 = and /x > 0. 

disappears completely. We have also observed additional minor effects that can play a role 
in the determination of the correct relic density. For example, a dependence of the higgsino 
component of the LSP on Vbl has been found. Again, this can lead to the disappearance 
of the Higgs funnel. Therefore, if SUSY is discovered with large tan this could again be 
interpreted as a lower limit on vr. 
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FIG. 13. Q.DA'lh'^ as a function of M1/2 for the choice of parameters niQ = 2.5 TeV, tan/3 = 30, 
Aq = and fi > 0. From left to right the different curves correspond to vbl = = [2 ■ 10^^, 
1.9 • 10^6, 1.8 • 10^6, 10^6, 5 • 10l^ 3 • 10^^ 2 • 10^^, 10^^] GeV. All curves reach their minimum, 
lower Ml/2 values would spoil electroweak symmetry breaking. The dashed horizontal lines show 
the values O^m/i^ = [0.1018,0.1228]. 

4- Focus point 

The focus point is also affected by the running of the parameters above the parity breaking 
scale. The observed dark matter relic density is obtained in this region thanks to annihilation 
into bosons, a process which is only effective when the higgsino component of the lightest 
neutralino is sufficiently large. This is provided by a small fi parameter, /i ~ Mi. 

As explained, MilmsusY) is typically smaller in the f2LR model than in the CMSSM. 
When one lowers the high energy VEVs Vbl and/or Vfi^, the resulting Mi{msusY) gets 
lowered as well, see eq. ( ITT]) . Therefore, the required tuning with the /i parameter can only 
be obtained by increasing Mi/2- A shift of the focus point region towards larger values of 
Ml/2 is thus expected. 

This expectation has been check numerically. Figure [T3l shows VtoMh"^ as a function of 
Ml/2 for the choice of parameters tjiq = 2.5 TeV, tan/3 = 30, Aq = and /i > 0. As vbl 
and are decreased, one needs to go to larger values of M1/2 in order to reproduce the 
observed relic abundance. We find that for Vbl = ^ij ~ 10^^ GeV one cannot make the 
annihilations sufficiently effective and the focus point disappears (for this particular choice 
of parameters). Note that in this figure all curves reach their minimum, lower M1/2 values 
would spoil electroweak symmetry breaking. We also would like to emphasize that the 
dependence on vbl and vb is very strong. This can be clearly seen from the non-negligible 
shift in Ml/2 obtained when one goes from vbl = vb = fncuT = 2 ■ 10^^ GeV to just 
= vr = 1.9- 10^^ GeV. 

We also found a clear dependence of QoAih'^ on the angles ^1^2, as shown in figure [TH 
This can be easily understood from the matching conditions at the SU(2)b breaking scale, 
eqs. ([6]) and ([7]). The running of the soft parameter fnjj^ is strongly affected by the quark 
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FIG. 14. Q^iMh'^ as a function of sin(02 — ^i) for the choice of parameters mo = 2.5 TeV, Mx/2 = 708 
GeV, tan /3 = 30, = 0, /i > and = vn = 10^^ GeV. The dashed horizontal lines show the 
values ^DMh^ = [0.1018,0.1228]. 



Yukawas Fg'^ at energies above v^i. This leads to different values for the /i parameter at 
the SUSY scale, affects the annihilation cross-section of the lightest neutralino and modifies 
its relic abundance. Another relevant parameter is rritop, which also has a strong impact on 
Tn?^^. In our runs we have used the value m^p = 171.1 GeV. However, we point out that the 
effect of taking a different value for m^p can be compensated by choosing proper values for 
the angles 61^2- 



B. Flavoured co-annihilation 



The VtLK model is in principle well suited for flavoured co-annihilation |33| . The tr can 
be made lighter by flavour contributions, making it possible to have co-annihilation with the 
lightest neutralino in points of parameter space where it would be impossible without flavour 
effects. Moreover, to obtain the correct relic density one must include flavour violating 
processes, like TrXi — 7/i or trTr rfi. 

It is possible to flnd regions in the mo-Mi/2 plane where Br(yU — > 67) and the other LEV 
decays respect the experimental limits while having flavoured co-annihilation. In addition 
to the usual flne-tuning that is required in the CMSSM to obtain the observed dark matter 
relic density, flavoured co-annihilation also requires to tune the neutrino mixing parameters 
in order to suppress Br(/i — )■ 67). In order to have large flavour effects in the slepton sector, 
to be able to reduce the mass of the tr and, at the same time, respect the experimental 
limits, one must flnd values for ^13 (the reactor angle) and 6 (the Dirac phase) that allow 
cancellations in Br(/i — )■ 67). We flnd that this cancellation is more effective if the mass of 
the lightest neutrino is non-zero. The r LEV decays are put under control by choosing rrio 
and Ml/2 sufficiently large. Figure [15] shows one example. On the left panel the cancellation 
of Br(/x — )■ 67) is obtained with 6 = n and 9is = 8°, while on the right panel with 612 = 9.5°. 
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FIG. 15. ^^-^o/i^ and Br(/i — t- 67) contour plots in the mQ-Mi/2 plane for 6 = tt, ^13 = 8^ (left) 
and = 9.5'^ (right). The other parameters are indicated in the figures, except for the fiLR 
parameters, vr = 10^^ GeV and vbl = 10^^ GeV. 



Note that this occurs for different values of mo and M1/2. 

Let us comment on some particularities of the f2LR model concerning cancellations in 
Br(/x — )■ 67). It turns out that, contrary to the minimal type-I seesaw model, where the 
cancellation in Br(/i — 67) occurs for a Dirac phase 5 = ir and for 6*13 ^ 5°, independently 
of the mo and M1/2 valueqj, for the f2LR model the cancellation occurs for different values 
of 6*13 depending on mo and M1/2. This fact might be puzzling, since the Yukawa structure 
that appears in the running of m| and m^ is the same as in the type-I seesaw. Therefore, 
if one wants to cancel this combination of Yukawas by choosing the right ^13, it should be 
possible to use the same value in type-I and f2LR, and this is not the case. 

The key point is that, in addition to m| and ml, one also has a soft trilinear term 
contribution, Tg, to Br(/i — )■ 67). This contribution is present in both amplitude coefficients, 
Al and Af>, and suppressed in the type-I seesaw due to its proportionality to charged 
lepton masses. However, in the QLK model this contributions turns out to be much larger. 
This can be understood from the RGEs. The 1-loop RGEs for Tg in the parity conserving 
74 1 contain terms of the type /V^l'^^ with k = 1,2. This will lead to two pieces 



regime 

after applying the matching conditions at Vf>, see eq. ([7]): the conventional f'^fYg, which is 
charged lepton mass suppressed, and the new f'^ fYy, which is not suppressed and gives rise 
to new contributions to Br(yU — )■ 67). 

This conclusion has been checked numerically, by studying the dependence of the differ- 
ent contributions to Br(/i — )■ 67) on ^13. We found that for large Aq the Tg contribution is 
dominant, which explains why Br(/i — )■ 67) does not cancel for 6*13 ~ 5° as in the minimal 
type-I seesaw. For = the Tg contribution is much smaller, and thus one finds an 

^ See reference j3l for examples. 
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approximate cancellation for 613 ~ 5°, as expected. 



V. LHC PHENOMENOLOGY 

In this section we investigate predictions for the LHC phenomenology of the flLK model. 
The effects of the heavy states are two-fold: they change the spectrum and induce flavour 
violating off-diagonal elements in the slepton and sneutrino mass matrices. The changes 
of the spectrum can lead to a potentially measurable difference between the states which 



are mainly selectron- and smuon-like |2ll . |75|-|77|. Here we are going to focus on the lepton 



flavour violation within the SUSY cascade decays as they occur for example in the decay 
chain 

Ql ^ qX2 qe-fi+Xi (22) 

As already stated above, in this model one has additional lepton flavour mixing for the 
R-sleptons as opposed to the usual seesaw mechanisms where the lepton flavour mixing is 
in the left sector only. 

In Figure [16] we show total rates for pp ^ X2 ^^'Xi with / 7^ /' as a function of M1/2 for 
^/s = 14 TeV, mo = 100 GeV, Aq = 0, tan/3 = 10 and /i > for different values of vbl = vr- 
Choosing Vbl = implies that one has approximately the same amount of flavour mixing 
for left- and right-sleptons. Here we have summed over all initial states containing squarks 
and gluinos and over all cascade decays leading to a X2- We have furthermore required that 
exactly one X2 decays lepton flavour violating and that no additional lepton occurs in any 
of the cascade decays of the corresponding event. For the calculation of the cross section 
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we have used the FASER-LHC package 78| which is based on the program PROSPINO 
Recent ATLAS 



53| and CMS j5J| data constrain already gluino and squark masses and are 
mainly interpreted in the context of the CMSSM. In our model the spectrum differs with 
respect to the CMSSM but we expect that for our choice of parameters gluino masses are 
excluded below about 1 TeV. For this reason we put dashed lines if the gluino mass is below 
this bound and full lines if the value is above. 

The neutrino Yukawa couplings are chosen such that Br(yU — 67) is close to its current 
bound of 2.4 ■ 10~^^. Here we assumed that the right-handed neutrinos are degenerate in 
mass. Let us first discuss the final states containing a r-lepton. We see an increase of the 
rates with increasing M1/2 until it reaches a maximum between 700 and 800 GeV depending 
on the value Vbl = and then drops with increasing M1/2. This behaviour is due to an 
interplay of three effects: (i) with increasing M1/2 larger values for the entries of are 
allowed as Br(/i — )■ 67) gets suppressed by the heavier spectrum. This implies larger flavour 
violating decay rates of X2. (ii) For fixed mo the left sleptons are heavier than X2 ^'^^ small 
Ml/2 and they have about the same mass for M1/2 in our examples between 500 and 600 
GeV depending on the value of vbl = vr- (iii) Heavier squark and gluino masses imply 
reduced cross sections. The shift of the maxima to higher values of M1/2 with decreasing 
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FIG. 16. o" X Br for — )• X2 ^ ~^ Xi ^ /U X (upper left plot), — )• ^ ~^ Xi ^ t X (upper right 
plot), pp ^ X2 ^ Xi T ^ (lower plot), as a function of M1/2 for ^/s = 14 TeV. The curves 
correspond to different values for vbl = vr: 10^^ GeV (blue), 5 x 10^^ GeV (red) and 10^^ GeV 
(green). The other CMSSM parameters have been fixed to: niQ = 100 GeV, Aq = GeV, tan/3 = 10 
and jj, > 0. The seesaw scale Ms has been adjusted at each point in order to get Br(// — ?• 67) close 
to but smaller than 2.4 • lO'^^^ j^g 

gions with a gluino mass below 1 TeV are denoted with dashed 

lines. 

Vbl = Vr is again a consequence of the modified spectrum: smaller values oi Vbl = Vr imply 
smaller gaugino masses and increased ratios of slepton masses over gaugino masses at the 
electroweak scale for fixed M1/2, which shifts the M1/2 value where neutralino decays into 
on-shell L-sleptons are kinematically allowed. For completeness we note that the total rates 
can go up to 10 (30) fb in case of the er (/ir) implying at most a few thousands events if 
the design luminosity of 100 fb"^ per year can be achieved. The reduced values for the er 
channel compared to the /ir is a consequence of neutrino physics as we require tri-bimaximal 
mixing. This leads to 'rn\^^^lm\ 2z — ''^e=,i3/''^e'=,23 ~ fT^soilf^atm where mgoi and iriatm are 
the solar and atmospheric neutrino mass scale, respectively. 

In case of the e-/i final state in Figure [16] we see a second maximum at lower values of 
Ml/2 in a region where the final states containing a r have rates which are two to three 
orders of magnitude smaller. The reason for this is a level crossing of the with Jx^ 
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FIG. 17. a X Br for pp ^ X2 ^ ^' ^ function of | sin(02 ~ ^i)|- The bands correspond to 
the fohowing /,/' combinations: e -fi (blue), e-r (green) and /i-r (red). The parameters have been 
taken as in Figure [71 A center of mass energy of 14 TeV has been assumed. 

leading to an enhancement of three body decays X2 ~^ e/xx? via virtual left sleptons. At 
the maxima we find that Br(x2 e"'"e~Xi) — Br(x2 — Br(x2 f^~^f^~Xi) — 

Br(x2 T^^'y^) l'^- These three body decays give the dominant contribution to the flavour 
violating signal. 

Finally let us comment on the impact of the parameters we have kept fixed so far: for 
moderate increases, mo shifts the maximum of the rates as the mass of the sleptons is 
increased. For sufficiently large mo two body decays into sleptons become kinematically 
forbidden implying negligible event rates. A variation of Aq has only a small impact as it 
leads mainly to a shift of the masses for the states which are mainly stau-like. Larger values 
of tan /3 reduce the lepton flavour violating signal as Br(yU — )■ 67) grows like tan^ /3 implying 
smaller allowed values for the flavour off-diagonal entries in the slepton mass matrices and, 
thus, reduced branching ratios for the lepton flavour violating neutralino decays. 

Up to now we have fixed the parameters such that Br(/i — )■ 67) is close to its experimental 
bound. In the following we depart from this by performing a scan varying in particular 
I sin(6'2 — 9i)\ because a strong dependence of the rare lepton decays on this quantity has 
been found, see Figure [71 In Figure [T71 we show the rates for pp x\ X ^' ^ ^ 

function of | sin(6'2 — 6*1)! for ^/s = 14 TeV and fixing the parameters as in Figure [71 The final 
states containing a r-lepton behave similar as the corresponding rare r decays. However, 
in case of the e-fi final state we find a lower limit of about 0.1 fb for this parameter set. 
The reason is, that above | sin(6'2 — 0i)\ — 0.3 the mass splitting between the selectrons and 
smuons get smaller the larger | sin(^2 ~ ^i)| is- This holds for both, left and right states. 
This leads to a decrease of Br(/i — 67) and at the same time to a slight increase of the 
lepton flavour violating signals at the LHC as can also be seen in Figure [181 where we show 
the LHC rate of the e/i final state as a function of Br(yU — )■ 67). 

In principle one might worry about the fact that the mass differences between the sleptons 
are partly smaller than the corresponding widths implying potential resonance effects (soi 
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Br(/i ^ e 7) 



FIG. 18. a X Br for pp — t- X — t- ^ e ii X for y/s = 14 TeV as a function of Br(/i — t- 67) 
The parameters have been taken as in Figure [71 The vertical hne shows the current MEG bound 
Br(/i 67) < 2.4 • 10-i2_ 

El|- We have checked that the principal features discussed above remain if one does not use 
the narrow width approximation, leading to the cascades x% ~^ ~^ ^^'Xi? but calculates 
the complete amplitudes using Breit-Wigner propagators for the sleptons. This shifts the 
results only very slightly. 



VI. CONCLUSIONS 

We have studied phenomenological aspects of a supersymmetric left-right model. Our 
discussion has centered on two aspects: (a) Changes in the low energy SUSY spectra with 
respect to CMSSM expectations and the resulting consequences for the relic density of the 
lightest neutralino, assumed to be the cold dark matter of the universe. And (b) Lep- 
ton flavour violation induced by the LFV entries in the neutrino Yukawa coupling matrix, 
required to explain the observed neutrino angles, and consequences for the LHC 

In the CMSSM^the lightest neutralino can have a relic density in agreement with the 
measured value 
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only in some very specific parts of parameter space. These well-known 
regions are (i) the stau co-annihilation region; (ii) stop co-annihilation region; (iii) the focus 
point line and (iv) the Higgs funnel. The modified running of the soft parameters in the f2LR 
model shifts the allowed regions. In general, the lower Vbl and are, the smaller these 
regions become until for certain values of these parameters (depending on the specific region) 
the DM allowed regions disappear completely. From an observation of a SUSY spectrum 
consistent with any of these regions one could infer lower limits on vbl and vr not far below 
the GUT scale for CMSSM boundary conditions. 

The model also allows for flavoured co-annihilation jssl. We have given some specific 
examples of parameter points in which it is possible to have large flavoured contributions 
to the co-annihilation cross section, despite the stringent upper limit on the decay /x — 67 
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recently published by the MEG collaboration 56 

We have shown that LFV decays of the Xo, can reach up to 20 fb~^ for the LHC at 
^/s = 14 TeV. All combinations of different lepton flavour final states could be large. These 
as well as the rare lepton decays show a strong dependence on the parameters of the model, 
noteworthy also on the high scale parameters vr and sin(^2 — (^i)- 
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Appendix A: Invariants 



Neglecting the contributions from the Yukawa couplings and working in the leading-log 
approximation, the soft mass parameters at the SUSY scale are computed to be 
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In equ. ( 1A2[) a sum over the index i is implied, with i = (Y, 2, 3) in the first and second 
terms and i = {BL, 2, 3) in the thircfl. The /* parameters are defined as 
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^ In order to avoid any possible confusion at /i = vr, we use the notation ai{vii) = {ay, 012, oi^){vR) for the 
couphngs right below vr and aiivn) = (asL, a2, a3)(ffl) for the couplings right above vr- 
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For the numerical studies we have fixed the GUT scale to be rriGUT = 2 ■ 10^^ GeV. The 
gauge couplings in the previous equations are determined by means of the general formula 



l-^6Mnf4 



(A6) 



using the experimental values at /i = mz as starting point. In addition, at yU = vr one must 
apply the following matching conditions 
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In the previous three equations, the gauge couplings on the right-hand side are the ones for 
the ^ < vr regime, whereas the gauge couplings on the left-hand side are the ones for the 
H> Vr regime. 

The beta coefficients used in the previous formulas are 
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Finally, the A{ and coefficients are 
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